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We summarise results on the infrared behaviour of Landau gauge QCD from the Green's functions approach 
and lattice calculations. Approximate, nonperturbative solutions for the ghost, gluon and quark propagators as 
well as first results for the quark-gluon vertex from a coupled set of Dyson- Schwinger equations are compared 
to quenched and unquenched lattice results. Almost quantitative agreement is found for all three propagators. 
Similar effects of unquenching are found in both approaches. The dynamically generated quark masses are close 
to 'phenomenological' values. First results for the quark-gluon vertex indicate a complex tensor structure of the 
non-perturbative quark-gluon interaction. 



1. Introduction 

The nonperturbative properties of the quark- 
gluon interaction are at the heart of the most 
interesting phenomena of QCD such as confine- 
ment and dynamical chiral symmetry breaking. 
On a fundamental level the details of the quark- 
gluon interaction are encoded in the one-particle 
irreducible quark-gluon vertex. Together with 
the quark and gluon propagators this vertex en- 
ters as a vital ingredient in model building and 
plays a key role in bridging the gap between the 
coloured fundamental degrees of freedom of the 
theory and the observed colour-neutral hadron 
states. On a phenomenological level the quark- 
gluon interaction manifests itself in the details 
of the quark- (anti-)quark potential. Relativistic 
effects play an important role in potential mod- 
els and the proper choice of the Lorentz struc- 
ture of the quark-antiquark interaction is vital to 
guarantee the agreement of the theoretical pre- 
dictions with the experimental data. Ultimately, 
a derivation of these structures from the under- 
lying quark-gluon interaction is mandatory. To 
this end one has to know the properties of the 
most fundamental Green's functions of QCD, the 



dressed propagator and vertex functions. 

The infrared behaviour of the propagators of 
Landau gauge QCD has been investigated exten- 
sively over the past years in both, lattice Monte 
Carlo simulations and the continuum Green's 
functions approach. Lattice simulations are the 
only ab initio method known so far and are by 
now precise enough to pin down these propaga- 
tors accurately in a large momentum range cen- 
tered around 1 GeV. In the deep infrared, how- 
ever, lattice results are inevitably plagued by fi- 
nite volume effects. In the continuum formulation 
of QCD the Dyson-Schwinger equations (DSEs) 
provide a tool complementary to lattice simula- 
tions. They can be solved analytically in the 
infrared. Furthermore numerical solutions over 
the whole momentum range are available by now. 
The truncation assumptions necessary to close 
the DSEs can be checked in the momentum re- 
gions where lattice results are available. In gen- 
eral results from DSEs have the potential to pro- 
vide a successful description of hadrons in terms 
of quarks and gluons, see |1I2I3| and references 
therein. 
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2. Landau gauge QCD 

The ghost, gluon and quark propagators, 
Dcijp), D^^{p) and S{p), in EucUdean momen- 
tum space can be genericahy written as 



S{p,li^) = 



(1) 
(2) 

(3) 



Here denotes the renormahsation scale and 
G{p'^,IJ?) and are the ghost and gluon 

dressing functions. The Dyson-Schwinger equa- 
tions for these dressing functions have been solved 
in their continuum formulation |4I5| as well as on 
a torus, i.e. employing periodic boundary condi- 
tions 15]. In the continuum they can be solved 
analytically in the infrared and one finds simple 
power laws. 



Z{p\^^^) 
Gip^n') 



(4) 



for the gluon and ghost dressing function with 
exponents related to each other. Hereby k is an 
irrational number, k = (93 - \/2101)/98 « 0.595 
which depends slightly on the truncation scheme 
|7I8 |. Note that k > 0, also shown in im- 
plies the dominance of the 'geometric' ghost de- 
grees of freedom in the infrared in agreement 
with the Kugo-Ojima confinement criterion and 
Zwanziger's horizon condition 10 l l.,12j . Re- 
cently these results from the Dyson-Schwinger 
approach have been confirmed independently in 
studies of the exact renormahsation group equa- 
tion P] . 

A combination of the ghost and gluon dressing 
functions can be used to define the nonperturba- 
tive running coupling ^ 



a{p^)=a{^,^)GHp^^^^)Zip',^,'). 



(5) 



No vertex function appears in this definition; a 
fact that can be traced back to the ultraviolet 
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Figure 1. The gluon dressing hmction Z{p'^) from 
DSEs ISEl and on the lattice ,17.18, . 



finiteness of the ghost-gluon vertex in Landau 
gauge. Note that the right hand side of this 
equation is a renormahsation group invariant, 
i.e. a{p^) does not depend on the renormah- 
sation point. A further renormahsation group 
invariant is the quark mass function, given by 
Af(p2) ^ B{p'^,fi'^)/A(p'^,^'^). From the power 
laws Q one finds that the coupling has a fixed 
point in the infrared given by 

^ i^r(3-2.)r(3 + .)r(i + .) ^ 
^ ' 6Nc r2(2-K)r(2K) 

for the gauge group SU(3). The infrared dom- 
inance of the ghosts can be used to show that 
a(0) depends only weakly on the dressing of the 
ghost-gluon vertex and not at all on other vertex 
functions [7j. 

The running coupling as it results from numeri- 
cal solutions for the gluon, ghost and quark prop- 
agators can be quite accurately fitted by the rel- 
atively simple function |14 | 



a(0) 



47r pVA^ 



qcD 



l+p2/A2 



QCD 



QCD 



Hp'/Hcu) pVA^cd - 1 



(6) 



with /3o = (lliVc - 2iV/)/3. Note that, fol- 
lowing ref. ^3], the Landau pole has been sub- 
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Figure 2. The gluon propagator Z{p'^)/p'^ from 
DSEs [SEl and on the lattice [TTIT^ . 
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Figure 3. The quenched and unquenched quark 
mass function M(p^) and the wave function 
Zq(p^) from the DSE approach compared to 
results from quenched lattice calculations |16| . 



tracted. The scale Aqcd is hereby determined by 
fixing the running coupling at a certain scale, e.g. 
as{Ml) = 0.118. 

The numerical solutions for the gluon dressing 
function Z(ji'^) in the continuum and on the torus 
are compared to the results of recent lattice sim- 
ulations .17,18, in fig. ^ The qualitative agree- 
ment with the lattice results is very good and the 
quantitative discrepancies in the momentum re- 
gion around 1 GeV well understood: gluon self 
interactions not contained in the DSE-truncation 
play an important role here. In the infrared both 
approaches agree nicely. An interesting qualita- 
tive difference, however, can be seen in fig. El 
where we show the propagator Z{p^)/p'^ instead 
of the dressing function. Both, the lattice results 
and the DSE-solution on the compact manifold 
tend towards a constant in the infrared, whereas 
the propagator from the DSEs in continuum for- 
mulation vanishes. Both DSE-solutions employ 
the same truncation scheme, thus the differences 
between them have to be attributed to the peri- 
odic boundary conditions and the finite volume 
on the torus. These properties of the compact 
space might also account for the constant gluon 
propagator found on the lattice. 



3. Dynamical mass generation 

The dynamical generation of quark masses can 
be studied in detail in the Dyson-Schwinger equa- 
tion for the quark propagator, Eq.Q. It is a 
genuinely non-perturbative phenomenon and re- 
quires a careful treatment of the quark-gluon in- 
teraction employed. In the DSE- framework of ref . 
jl4| an ansatz for the quark-gluon vertex has been 
chosen that satisfies important constraints: it is 
genuinely nonperturbative in nature, guarantees 
multiplicative renormalisability in the quark DSE 
and has the correct limit in the perturbative mo- 
mentum domain. Its structure is such that it fac- 
torizes in an Abelian and a non-Abelian part, 

r.(g, k) = q, k) W-'""^\p, q, k), (7) 

with p and q denoting the quark momenta and 
k the gluon momentum. For the Abelian part 
Vj/*"^' carrying all the tensor structure the Curtis- 
Pennington construction has been chosen 
which has been used frequently in QED. This 
choice has the advantage that vector as well as 
scalar tensor components are taken into account 
and their effect on the quark propagator can be 
studied. In fact, as detailed in refs. [20], the 
inclusion of the scalar parts of the vertex is ca- 



4 



2.5 
2 

nI.5 



0.5 



1 


' □ N,=b, Bowman et al. 




Nj=2+1 , Bowman et al. 




i — DSE, N,=0 




% -- DSE, N,=3 (chiral) 




... DSE, N,=2+1 (u/d,s) - 












.1.1.1.1 



"0 1234567 

P [GeV] 

Figure 4. The quenched and unquenched gluon 
dressing function from the DSE approach ^1] 
compared to results from unquenched lattice cal- 
culations [TH]. 



pable of changing the analytic structure of the 
quark propagator dramatically, leading to a pos- 
itive definite spectral function for the quarks. 

The non-Abelian part W^'^''^\p,q,k) of the 
quark-gluon vertex Q contains ghost-dressing 
factors as (partly) implied by its Slavnov- Taylor 
identity. The combined dressing of the quark- 
gluon vertex and the full gluon propagator in 
the quark DSE then provides enough interaction 
strength to generate dynamical quark masses of 
the order of M{0) w 350 MeV in agreement with 
phenomenology. Furthermore a chiral conden- 
sate of the order (-(vi/*))-i/3 ^ 300 MeV has 
been obtained, which strongly favours the stan- 
dard counting rules in chiral perturbation theory. 

The results for the quenched quark mass func- 
tion M{p^) and the wave function Zq{p^) are 
compared to the quenched lattice results of refs. 

in fig. 121 The overall qualitative and quan- 
titative agreement between both approaches is 
very good. The DSE results are within the 
bounds given by the two different formulations 
of fermions on the lattice. 

Including the backreaction of the quark- 
propagator on the ghost and gluon system leads 
to a coupled set of three Dyson-Schwinger equa- 



tions for the propagators of QCD. These equa- 
tions have been solved in jj^ and allowed a pre- 
diction of possible effects of unquenching QCD 
on the propagators. As can be seen from fig. |2| 
including Nf = 3 chiral quarks in the gluon DSE 
hardly changes the results for the quark propaga- 
tor. The slight enhancement in the perturbative 
tail of the mass function is in perfect accordance 
with the expected change of the anomalous di- 
mension. The chiral condensate is nearly unaf- 
fected. It will be interesting to compare these 
results to unquenched lattice calculations when 
available. 

Unquenched lattice results for the gluon prop- 
agator including the effects of two light (up-) and 
one heavy (strange-) quark have been published 
recently |18j and are compared to the correspond- 
ing results from our DSE-approach in fig. ^ The 
screening effect from the quark loop is clearly vis- 
ible in the lattice results for momenta p larger 
than p = 0.5 GeV: the gluonic self interaction be- 
comes less important in this region and the gluon 
dressing increases. Although to a somewhat less 
extent, this effect can also be seen in the DSE- 
approach. As mentioned already above not all 
effects from the gluonic self interaction are ac- 
counted for in the DSE truncation. When this 
part of the gluon interaction becomes less domi- 
nant in the unquenched case, both the lattice and 
the DSE-approach agree very well on a quantita- 
tive level, provided similar bare quark masses are 
taken into account. In the chiral limit the screen- 
ing effect of the quark loop becomes even stronger 
as can be seen from the DSE-results in fig. 0] 
This is certainly expected as the energy needed 
to create a quark pair out of the vacuum becomes 
smaller with decreasing bare quark mass. 

Both, the lattice calculations and the Green's 
functions approach agree in the fact that un- 
quenching does not affect the extreme infrared of 
the ghost and gluon propagators. Again, this is 
easily explained from dynamical chiral symmetry 
breaking: there is not enough energy to generate 
a quark pair from the vacuum below a certain 
threshold, and the quark degrees of freedom de- 
couple from the Yang-Mills sector of the theory. 



5 



chiral quarks; (p,)'=p , (p^) =(2p) , (p,-p,) =(3p) 



2 2.5 7 




1 1.5 

P (GeV) 



2.5 r 



m = lis MeV; (pj)-=(pj)-=:(p) 



: 5 



A,| (lattice) 
(4p^) (lattice) 
(-2p)Xj (lattice) 

(DSE) 
(4p') (DSE) 
(-2p) X (DSE) 




1 2 

P (GeV) 



Figure 5. The lea(iing tensor structures of the 
quark gluon vertex from the vertex DSE with chi- 
ral quarks. 

4. The quark-gluon vertex 

As mentioned in the introduction, the quark- 
gluon vertex plays an important role in both, 
understanding confinement from the quark-gluon 
interaction and providing the bridge between 
coloured quarks and giuons, and their colourless 
bound states, the hadrons. Therefore it is de- 
sirable to determine its properties in a selfcon- 
sistent calculation. As a first step in this direc- 
tion we present a calculation of the components 
of the quark-gluon vertex from its DSE (also see 
Ref.|21|). Our calculation is not selfconsistent 
yet: we use the quark, ghost and gluon propa- 
gators presented in the last section as input on 
the right hand side of the vertex DSE without 
backcoupling the vertex to the quark equation. 
Furthermore we only employ the 7^-part of the 
vertex in all loops of the vertex DSE. Neverthe- 
less we believe that our results are meaningful as 
the propagators used are close to their respective 
counterparts on the lattice. Furthermore from 
the structure of the vertex DSE one finds indi- 
cations that selfconsistency effects could play a 
major role in the infrared but are small in the 
medium and large momentum regime. 

With these caveats in mind we present our re- 
sults in figs. [SI and IHI The complete vertex can be 
decomposed in a basis of twelve tensor structures, 

- f (Eti M + Ell nT, m) EH- Here 



Figure 6. Tensor structures of the quark gluon 
vertex from the vertex DSE with a massive quark 
compared to lattice calculations )22j . 

we show only five of these components, namely 
the ones multiplying the tensors 

Li ^ = 7m 

L2 = -{i>i+ h){pi+P2)tJ. 

T5 f, = icrf,i.{P2 - PiT 

M = -l^L(J^^PlP2~ i>iP2 hPi M 

with the quark momenta pi and P2- The (di- 
mensionless) four leading structures in the chiral 
limit are given in fig.|Sl Although the Li- piece is 
dominating, sizeable admixtures from other ten- 
sor structures occur. This is also found in a re- 
cent model study |22I and it is expected, as the 
quark-gluon vertex provides one of the underly- 
ing structures from which the rich structure of 
phenomenological quark potentials (see e.g. |24j) 
should be ultimately derived. 

Our results for a massive quark of mo = 115 
MeV are compared with the available lattice re- 
sults in fig.ini Even without selfconsistency in the 
DSEs we obtain very nice agreement for the vec- 
tor and the scalar tensor components, Li and L3. 
On the other hand, the infrared behaviour of the 
L2-piece is in obvious disagreement. Firm con- 
clusions, however, cannot be drawn before lattice 
data with smaller errors as well as a selfconsistent 
DSE solution are available. 
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5. Outlook 

In the last years a consistent picture of the 
infrared behaviour of Landau gauge QCD has 
emerged. Evidence from the Green's functions 
approach and lattice calculations suggests that 
the Faddeev-Popov determinant (i.e. the ghost 
degrees of freedom) is dominant in the infrared 
and provide the long range interaction in QCD. 
The gluon propagator is finite or even vanishing 
in the infrared. Taken at face value this behaviour 
provides a problem for confinement: Simple one- 
gluon exchange with a bare quark-gluon vertex 
cannot account for the linear rising potential be- 
tween two static colour sources. Infrared singu- 
larities present in the dressed quark-gluon vertex 
could resolve this issue. Whether this is indeed 
the case or a more sophisticated mechanism is at 
work remains an open problem for the future. 
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